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ABSTRACT 

We discuss the latest CCFR determination of the strange sea density of the proton. 
We comment on the differences with a previous, leading-order, result and point out the 
relevance of quark mass effects and current non-conservation effects. By taking them 
into account it is possible to solve the residual discrepancy with another determination 
of the strange quark distribution. Two important sources of uncertainties are also 
analyzed. 
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1 Introduction 



Until few years ago the problem of the strange sea distribution was rather controversial, 
despite the impressive amount of knowledge on the internal structure of nucleons ac- 
cumulated in the last decade. The best available fits in 1993 [|l], differed by almost a 
factor 2 for the strange density. As we pointed out in @, f§> this puzzling situation was 
mainly determined by an erroneous interpretation of the experimental results, which 
ignored important physical effects already investigated in || [|. 

There are two ways to extract the strange sea distribution s(x). The first method 
consists in subtracting the structure functions F 2 measured in muon (/i) and neutrino 
(u) deep inelastic scattering (DIS). We call this the v — \i determination of s(x). The 
second, more direct, method consists in selecting iDIS events with charm production: 
the signature of these events is the presence of opposite-sign dimuons in the final state. 

When the NMC //DIS data and the CCFR i/DIS data g | made both these 
determinations possible, it was found, rather surprisingly, that the two results for s(x) 
were largely different (see however the anticipations in || [|): the strange density 
extracted from dimuon data was considerably smaller than the one obtained by sub- 
tracting z/DIS and /iDIS data. 

In trying to solve this discrepancy all the available global parton parametrizations 
ran into serious difficulties. The CTEQ group, trusting the v — fi result, was led to 
advocate a very large strange sea content k = 2S/(U + D) ~ 0.9 (S = J dxxs(x), 
etc.), i.e. an almost 5*?7(3)-symmetric light sea. At the same time, however, the 
CTEQ strange distribution lied well above the dimuon data. On the other hand, the 
MRS D' Q fit stuck to the more conservative value k = 0.5 (which was an input) but still 
overshot the dimuon data, while lying below the v — \x result. (We shall see a 'posteriori 
that the MRS compromise between the two data sets and the value of k ~ 0.5 are 
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-accidentally- much closer to the truth than the CTEQ picture.) 

The solution to the strange sea puzzle proposed in ||, |4[] was very simple: there is no 
real puzzle because the two determinations of s(x) measure in fact different quantities, 
none of which coincides with the true strange density. This is due to the fact that, 
up to moderately large Q 2 (of order of 30GeV 2 or so, i.e. not much above the charm 
threshold), the relevant diagrams for charm production are the vector-boson-gluon 



fusion subprocesses ]nj. These are, in the common massless QCD terminology, next- 
to-leading order diagrams and hence they are often mistakenly forgotten as if they were 
subleading corrections. When the gluon-fusion diagrams are taken into account, two 
effects arise which explain why the two determinations of s(x) do not really provide 
s(x) (at least directly, as a naive leading-order analysis would suggest) and should 
indeed give different results. They are: i) quark mass corrections; ii) non conservation 
of weak currents (yielding large longitudinal contributions). These effects are calculable 
in perturbative QCD |5|, [| and are relevant up to moderate Q 2 values (not much larger 
than the heavy quark mass scale). We called them non-universality effects because 
they make the heavy flavour contributions to iaDIS and to /iDIS structure functions 
intrinsically different. 

Near threshold, massless QCD is inappropriate to describe heavy quark production 
and mass contributions must be kept. The gluon-fusion process is the dominant one. 
Since the mass thresholds in the transitions W + s — > c, 7* — ► ss(cc), Z° — > ss(cc), are 
different, and since the longitudinal structure functions in weak DIS are larger than in 
electromagnetic DIS because weak currents are non conserved, we expect |5|, |j] 

where we denote by the sc contribution to the z/DIS structure function F2 with 
the electroweak coupling factored out (differently stated, F 2 I/,SC reduces to x(s + c) at 



leading order), and we use an analogous notation for the other quantities. 

Due to the importance of the gluon-fusion processes, we must expect a considerable 
difference between a leading order (LO) and a next-to-leading (NLO) extraction of 
s(x) from the dimuon experiment 0, |J. This expectation proves correct. The CCFR 



collaboration has recently released [i 1] a new determination of s(x) based on an analysis 
of the dimuon data which takes into account the gluon-fusion processes. The 'new' 



strange density [1 Ij is considerably larger than the 'old' one and partially bridges 
up the gap with the v — /i result. 

The purpose of this paper is to discuss and clarify our present knowledge of the 
strange density. We shall comment on the recent CCFR determination of s(x), on its 
difference with the previous one, and on the important physics behind such a difference. 
On the quantitative side, we shall show that the residual discrepancy existing at small 
x between the CCFR NLO strange density and the v — fi data is easily accounted for 
by the non-universality (or, differently said, by next-to-leading order) effects related to 
the gluon fusion processes. Practically, our calculations endorse an MRS-like strange 
sea density |T2| with a value k ~ 0.5. We shall also discuss some subtleties and 
some sources of uncertainties in the analysis of the neutrino data, and propose a more 
convenient way to present the dimuon results. 



2 The strange density from dimuon data 

The most direct mechanism to measure the strange density is charm production in 
charged current neutrino deep inelastic scattering. 
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The z/DIS cross section reads |13 
d 2 a vN G 2 m N E u 



dxdy 



7T(1 + 



Mr, 



xy'Fr (x) + (l- y)E»» (x) + (y - *-)xF^ (x) 



(2) 



If we restrict ourselves to charm excitation, only the transitions W + d — > c, W + s — > c 
contribute to the structure functions. At leading order we have to consider only the 
quark scattering terms [] and the cross section for charm production in z/DIS can be 
expressed in terms of the LO parton densities as 

d 2 a(uN -> cX) G 2 m N E v 



dxdy 



Trfl 



x{[u(x,Q 2 )+d(x,Q 2 )] \V cd \ 2 + 2s(x,Q 2 )\V cs \ 2 } (3) 



where V c d, V cs are the Cabibbo-Kobayashi-Maskawa matrix elements. Since V cs ^> V cc i 
the measure of a(uN — > cX) provides an excellent determination of xs(x). 

In order to take into account effects connected with the non negligible mass of 
the charm quark, it has been customary for many years to adopt the slow-rescaling 
procedure (this is what the CCFR collaboration also did in |§). 

In practice, the slow rescaling method consists in replacing Bjorken's x by the new 
variable £ = x(l + m 2 /Q 2 ), which is (naively) expected to take into account the effects 
of the heavy quark mass. If the Callan-Gross relation is enforced in terms of £, namely 
F 2 (£) = 2£Fi, an extra factor appears in the z/DIS cross section which then reads 

d 2 a(uN -> cX) G 2 m N E u 



d^dy 



tt(1 + 



01 
Mr, 



X 



£ { Q 2 ) + Q 2 )] \V cd \ 2 + 2 S (C, Q 2 )\V CS \ 2 } 1 1 + 



2m N E v £ 

Slow rescaling lacks a solid theoretical foundation. It is a sensible method if one 
considers only the quark scattering term. In this case, in the W + s — > c transition, the 
s quark is taken on shell, as usual in the parton model, and its mass is neglected with 
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2 For the sake of clarity we stick hereafter to the terminology of Ref. £3, which is used by the 
CCFR collaboration in their analysis. 



respect to the charm mass: if we call £ the fraction of proton's momentum carried by 
tha strange quark, the substitution x — > £ follows straightforwardly. 

However, considering higher order, gluon fusion, diagrams (which are actually the 
dominant ones) it is no longer possible to assume that the s quark is on shell: the 
whole procedure thus breaks up and a more sophisticated treatment is called for. Of 
course, one could still think that slow rescaling mimicks rather faithfully the real world 
and that it accounts for quark mass effects in an effective way. This attitude has been 
quite popular and it is still widely believed that slow rescaling is at least a very good 
approximation. One of the conclusions of the first, leading order, CCFR analysis was 
that the data supported the slow rescaling model of charm production (although with a 
very small charm mass, m c ~ 1.31 GeV 2 ). The NLO analysis has completely reversed 
the situation: a large difference is found between the LO determination with slow 
rescaling and the new one, which is on a better theoretical ground and leads to a more 
realistic charm mass, m c ~ 1.70 GeV 2 . The first conclusion we can draw from the new 
results is indeed that the slow rescaling, which is intrinsically a LO procedure, fails to 
provide a realistic picture of charm production and therefore must simply be abandoned 
(this criticism was anticipated in ||, well before both CCFR determinations). 

In the region of small and moderate values of Q 2 (where most of the CCFR data 
lie), i.e. not much above the charm threshold m 2 ,, it is not legitimate to retain only 
the quark scattering LO diagrams. As a matter of fact, near threshold, the whole 
contribution of a heavy quark to structure functions^ is given by the vector-boson- 
gluon fusion process (Fig. 1), which are conventionally classified as a next-to-leading 
order term (although the leading term is not the dominant one). 

In the case at hand we have for F 2 , in a formal notation (® means convolution) 

3 Remember that, in charged current ^DIS, probing the s quark means at the same time exciting 
the c quark. 
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01 

F? = (^)g®C 2 (W + g->sc), (5) 

where g(x, Q 2 ) is the gluon distribution and C 2 is the unsubtracted Wilson coefficient, 
i. e. the full cross section for the W + g — > sc process, which is made of the two insepa- 
rable diagrams shown in Fig. 1. 

It is worthwhile to draw a comparison with another approach [ 14 1 . In the formalism 



of Ref. ]l4| the sc contribution to structure functions is expressed as (we omit all the 
electroweak couplings) 

Ff = x(s + c) - CT htT + (^)g® C 2 (W + g - sc) , (6) 

where s(x, Q 2 ) and c(x, Q 2 ) are NLO parton densities, and C| ubtr is a subtraction term 
proportional to a s \og{^ 2 /m 2 ) (fif is the factorization scale). In this approach the 
logarithmic term in the Wilson coefficient, which explodes at Q 2 much larger than m 2 , 
giving eventually rise to collinear singularities, is subtracted out; at the same time, 
the quark scattering term is introduced so that, when the physical scale is very large 
compared to the heavy quark mass, the massless QCD parton model is regained. 

Up to Q 2 of order 10 m 2 the two formulas for F| c given above, eqs. (|5|) and (|GJ), are 
equivalent ||14|| . In fact, the two extra-terms in eq. @ are approximately equal and 
cancel out. All the relevant physics (mass effects and current non conservation effects 
producing large longitudinal structure functions) is thus contained in the W-gluon 
fusion diagrams. 

When these are taken into account the relation between the i^DIS cross section 
with charm production and the parton densities is much more involved than eqs. (0) 
or (U). In fact, the simplicity of the leading order formulas for d 2 a{vN — > cX)/d£d?/ is 
determined by the (fortuitously) virtuous combination of two elements: i) the relation 
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(valid for the structure function components relevant to charm excitation) 

F£ N = 2xF? N = xF» N = x{u + d + 2s) , (7) 

(or the corresponding one in terms of £), and ii) the consequent, accidental, cancellation 
of all factors in y in front of the structure functions. None of these two circumstances 
occur in the next-to-leading order case (we keep using this terminology although we 
stressed above that the Wg fusion contribution is the dominant one and cannot be, 
strictly speaking, called next-to-leading). The expression on the r.h.s. of eq. (Q) 
is therefore much more complicated and, of course, the slow rescaling substitution 
becomes completely meaningless. 

Hence it is expected, from a theoretical viewpoint, that there should be a large 
difference between the leading order strange density determined from @ and the next- 
to-leading order strange density extracted from (Q) upon use of eq. (||) . This is indeed 
what has been recently found |Tl| in the new, NLO, CCFR analysis of the dimuon 



data. As one can see in Fig. 2, at moderate Q 2 values, the difference between the LO 
result for xs(x) and the NLO result is rather large (more than 50% at Q 2 ~ 10 GeV 2 ). 
This is clearly not a mere higher-order correction: the LO analysis hides and neglects 
all the important physics contained in the H^-gluon fusion diagrams. The slow rescal- 
ing method, while theoretically ill-founded, is not even an effectively successful way 
to account for heavy quark masses. Moreover, mass effects are not the whole story. 
Effects of current non conservation in z/DIS scattering are quantitatively as much (or 
even more) important |§. In z/DIS the ratio R = o^jaT in the heavy quark sector 
near threshold is larger than 1. Moreover it is different in charged current processes 
(where both the vector and the axial currents are not conserved) and in neutral cur- 
rent processes (where only the axial current is not conserved). The large longitudinal 
contribution leads to a strong violation of the Callan-Gross relation. Neither the use 
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of this relation nor its correction by a value of R taken from electromagnetic scattering 
(as it was done in the LO analysis ||) are then legitimate. The inclusion of gluon 
fusion diagrams allows considering all important physical effects automatically and in 
a QCD computable way. 

The new CCFR strange quark distribution is in much better agreement with the 
MRS-A fit (see Fig. 2) and leads to a strange sea content re = 0.477 ± 0.050 
at the average (Q 2 ) ~ 22GeV 2 . There is a (10 — 15)% error on the result due to 
the factorization scale uncertainty (see sec. 4). The measured quark mass is m c = 
1.70 ± 0.19 GeV 2 and is larger than the unrealistic value m c = 1.31 GeV 2 found in the 
LO analysis of Ref. || because of the slow rescaling method adopted there. 

Finally, notice that, in passing from LO to NLO, a variation is expected (and 
detected ||) also for the non strange distributions, which are Cabibbo suppressed in 
the cross section ||, because the u and d contributions to structure functions are mixed 
with the charm contribution in the W + g — > uc, dc processes. 



3 Comparison with another determination of the 
strange density 

There is another way to extract the strange sea density from deep inelastic scattering. 
It combines measurements of z/DIS and /iDIS structure functions: we shall call it the 
v — fi determination of s(x). 

Let us resort once more to the parton model or, equivalently, to leading order QCD. 
The decompositions of F 2 for muonic and neutrino probes and an isoscalar nucleon are 

F£ N = —x{u + u + d + d) + ^rxc + ^-xs , (8) 
lo 9 9 
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F% N = x{u + u + d + d + 2s + 2c). (9) 

In eq. (§) we used s = s, c = c and eq. @ refers to charged current scattering. 
By appropriately combining F% N and F£ one can select the (leading order) strange 
distribution, under the assumption c<Cs: 

\f» n (x, Q 2 ) - 3F£ N (x, Q 2 ) = xs(x, Q 2 ) . (10) 

D 

It is immediately evident what is the main experimental difficulty with this deter- 
mination. The quantity on the l.h.s. of eq.(|H]|) is obtained by subtracting data from 
two different experiments and is very sensitive to the relative normalization. Besides, 
data are not taken at the same x and Q 2 values. Large uncertainties thus arise in the 
v — fi difference. 

On the other hand, the v — \i determination presents at least two advantages. First 
of all, on the theoretical side, its parton density content is simpler and much easier to 
reconstruct than that of the quantity in curly brackets in eq. @. Second, when F% N 
is measured there is no spurious, acceptance-dependent, separation of the two gluon 
fusion diagrams in Fig. 1, as in the dimuon measurement (see next section), and the 
whole sc contribution to the structure function is determined. 

Now, after the first (LO) CCFR determination of the strange density came out || it 
was clear that there existed a big discrepancy (see Fig. 2) with the result from the v — \i 
difference, obtained by combining NMC [7j and CCFR || data, although the latter 
was affected by large errors. As we have already recalled, this discrepancy was mainly 
due to the leading order analysis of the dimuon data which hided important physical 
effects. With the new CCFR determination of xs(x) which includes the contribution of 
the gluon fusion diagrams the situation is considerably less puzzling. Still, there seems 
to be a residual discrepancy with the v = fi result, which deserves some explanation 
(see Fig. 2). This is simple if one remembers that the non-universality effects taken 
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into account in the NLO dimuon result are not included in the LO formula (|T0| ) on 
which the identification of the difference (5/6)^^ — 3 with the strange density 



is based. Otherwise stated, NLO effects invalidate eq. (|T0f) which must be replaced by 



^FZ N -3F£ N = ^F| c -~Ff , (11) 
where (we are interested in the Q 2 region around 10 GeV 2 ) 

F? = (-^)g®C 2 (W + g^sc) } (12) 

Ff s = x(s + s) + (-^)g®Cr ssless (l*g-*ss)^x(s + s) . (13) 

Z7T 



The second equality in eq. flT3"|) is valid when the physical scale Q 2 is sufficiently higher 
than the strange threshold. 

Were the non-universality ratio F^/F^ equal to 1/2, as in the LO case, eq. (|Tl|) 
would reduce to eq. (|10"1) . However, for not too large Q 2 , this ratio is largely different 
from 1/2 (see || |J and, for a more systematic study, fl5|| ). Therefore the difference 
(5/6) F% N — 3F 2 MAr does not coincide with the strange quark distribution xs{x). We 
evaluated the r.h.s. of eq. (|TT|) by resorting to eqs. (|12||i"3|) , with factorization scale 
fi 2 = m 2 , and using the MRS-A parton densities [[Tj|. Since the MRS-A strange 
density reproduces rather well the CCFR NLO data, our calculation will also clarify 
whether there is a real contradiction between the dimuon and the v — \x determinations. 
The result for the v — fi difference ([11]) is displayed in Fig. 3 and compared to the 
data. The good agreement found shows that the dimuon and the v — \x determinations 
are compatible in the whole x range, and is a check of the goodness of the MRS-A 
parametrization of xs(x). Notice that, accounting for NLO effects, the v — fi difference 
turns out to be larger than xs(x) (dashed line) at low x and smaller at high x, with a 
crossover at x ~ 0.07. 

With the (important) caveat illustrated in the next section, we have now a trustwor- 
thy picture of the strange density. If we believe in the recent CCFR determination and 
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we assume the reliability of the v — /j, data, no discrepancy whatsoever is at present 
detected and all experimental determinations converge to an s(x) well fitted by the 
MRS curve with a conservative value of k (~ 0.5). 



4 Uncertainties in the extraction of the strange den- 
sity 

We have seen that the strange density recently determined by the CCFR Collaboration 
has been obtained by a next-to-leading order QCD analysis of the dimuon data, which 
considers all the relevant effects previously overlooked. The result, however, is affected 
by a relatively large inherent uncertainty (more than 10%). 

In this section we want to discuss two important sources of systematic uncertainties 
in the extraction of the strange sea density: one is specific of the dimuon determination, 
the other is more generally related to the kinematical region considered, close to a heavy 
quark threshold. 

The first correction was already discussed in 01], where it was pointed out that 
in the first CCFR extraction of the strange density an important acceptance effect 
had been neglected, namely the experimentally different weight of the two diagrams in 
Fig. 1 due to the energy cut on the second muon. 

In order to limit the background (mostly due to kaon and pion decays) a lower 
cutoff is set on the momentum p M2 of the muon coming from the semileptonic decay of 
charm (p^ 2 > 5 GeV) . This reflects itself into a cut on the momentum of the produced c 
quark (we shall call z the fraction of the light-cone momentum of the W boson carried 
by the charm). Given that the low-z region is dominated by the subprocess where 
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the gluon splits into a cc pair with the s produced in the ^-hemisphere (u channel), 
and, viceversa, the high-z region is dominated by the process where the gluon splits 
into an ss pair with the c produced in the ly-hemisphere, it is clear that the cut on z 
introduces an acceptance correction on the final result for the strange density. A way to 
compute such correction is to look at the ^-distribution aw g {.z) for the W^-gluon fusion 
process @]. This is strongly asymmetric: the two peaks at z — > (backward peak) 
and at z — > 1 (forward peak) are not delta-like at small Q 2 (they tend to become 
so only at asymptotically large Q 2 ) and have a different Q 2 evolution: the forward 
peak appears sooner. The p M2 cut can be taken to generate an effective cutoff z c on 
z. The jy-gluon fusion contribution to the structure function F£ c really measured 
is then proportional to J z awg(z)dz and hence smaller than the whole contribution 
Jq 1 aw g (z)dz. Were z c known, the correction would be theoretically predictable. On 
the other hand a possible uncertainty on z c may be a non irrelevant source of error on 
the extracted strange structure function. In Fig. 4 we show the results of the model 
of Ref. Q for F% c with two choices for z c : z c = 0.5 (solid curve) and z c = 0.8 (dashed 
curve). The area below these two curves amounts to about 75% and 55%, respectively, 
of the whole integral (corresponding to z c = 0). 

An experimental evaluation of this acceptance correction can be performed by fold- 
ing the z-cross section (Twg(z) with an empirically determined acceptance function, 
which is zero at z —>■ and rises to 1 at large z (the method sketched in the previous 
paragraph corresponds to a step-function choice for the acceptance curve). In its lat- 
est analysis O] , the CCFR collaboration has carried out such computation and found 
that the acceptance correction is 60% with an estimated error of 10%. These values 
correspond approximately, in our approach, to the situation depicted in Fig. 4 and 
confirm the importance of the effect (and also, incidentally, the educated guess on the 



13 



acceptance uncertainty made in Q). 

A more relevant (and fundamental) uncertainty on the strange quark density comes 
from the fact that at Q 2 not much below a heavy quark threshold it is generally unsafe 
to extract the distributions of the vector-boson-gluon fusion products, instead of their 
contributions to structure functions. 

To explain this point, let us look for instance at the expression (|]) for F| c . Since 
the jy-quark fusion terms are negligible with respect to the W-gluon fusion ones, 
the (anti) quark distributions are contained only in the quark scattering term. This 
is small near threshold and undergoes a subtle and rather precise cancellation with 
the subtraction term: therefore its extraction is a delicate matter. More importantly, 
the quark scattering contribution is strongly dependent on the factorization scale, as 
it is theoretically predicted [H and experimentally observed [IT | : most of the overall 
(10 — 15)% uncertainty in the NLO CCFR result for xs(x) comes from the arbitrariness 
in the choice of /i 2 . By contrast, the structure function F| c is rather stable against 
various choices of the factorization scale and is therefore the best quantity to explore, 
at least as far as the quark scattering terms do not become the dominant contribution 
- which happens at Q 2 well above threshold. 

It would thus be desirable, at least near heavy quark thresholds, to get from exper- 
iments data on structure functions instead of data on parton distributions. 



5 Conclusions 

Let us summarize the main points of this work. 

We have by now a much better knowledge of the strange content of the nucleon, 
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coming mostly from neutrino deep inelastic scattering. 

The next-to-leading order determination of xs(x) performed by CCFR supersedes 
the leading-order one: the latter should be recorded as a result which has little to do 
with the real world. It relied on two assumptions: i) the gluon fusion diagrams are neg- 
ligible corrections; it) the quark mass effects are accounted for by slow rescaling. In the 
Q 2 region of present experimental interest, we pointed out that: the former assump- 
tion is simply wrong (as the explicit calculations show); the latter is ill-established 
and the slow rescaling procedure does not even mimick the correct treatment of heavy 
quarks. Not far from threshold, the next-to-leading order term is not a correction of 
the leading-order term and one should be very careful in using such a terminology 
which could induce into dangerous misunderstandings. 

The discrepancy between the dimuon and the v — /i results for xs(x), which was 
worrisome at the time of the first CCFR determination, is now solved. The higher- 
order analysis correctly takes into account the effects which were the physical cause of 
such an apparent puzzle (we dubbed them "non-universality effects"): different mass 
thresholds in neutral and charged current DIS and large longitudinal contributions to 
z/DIS structure functions. In other terms, the dimuon data and the v — \x data measure 
different quantities, which coincide only when the two classes of effects just mentioned 
are neglected, i.e. only when a LO analysis is performed. By accounting for these 
non-universality (or NLO) effects [i.e. for the non negligible difference between 2F| C 
and F 2 CC + -F 2 SS ), we have explicitly shown that even the residual gap with the new 
dimuon data at small x is fictitious. 

Although greatly improved, our present knowledge of the strange distribution is 
not yet free from uncertainties. The NLO extraction of the strange sea density from 
dimuon data near charm threshold is intrinsically unsafe, because the quark scattering 
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term (that is, the parton density) is a small contribution subjected to cancellation by 
the subtraction term and, at the same time, has a relatively strong dependence on 
the factorization scale. This dependence weakens if one considers the whole structure 
function (including the dominant gluon fusion term). Thus, at least at moderate values 
of Q 2 , the strange and charm structure functions (and not their parton distributions) 
should be experimentally extracted to be object of theoretical study. 
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Figure captions 



Fig. 1 P^-gluon fusion process for charm production: a) u channel diagram; b) t channel 
diagram. 

Fig. 2 The strange quark distribution xs(x). The boxes are the LO CCFR determination 
U at Q 2 ~ lOGeV 2 . The circles are the quantity (5/6)^ - 3 F£ N at Q 2 ~ 
7GeV 2 , obtained by combining NMC and CCFR data [§. The shaded area 
represents the new NLO CCFR result [O] at Q 2 = 10 GeV 2 . The continuous line 



is the MRS-A fit for xs(x) at Q 2 = 7 GeV 2 



Fig. 3 The difference (5/6) F% N — 3 F% N at Q 2 ~ 7Ge\^ 2 . The circles are the experimen- 
tal results obtained by combining NMC and CCFR data ||. The solid line 
is the next-to-leading order QCD prediction described in the text. The dashed 



line is the MRS-A fit pj] for xs(x) at Q 2 = 7 GeV 2 



Fig. 4 The charm-strange structure function F| c at Q 2 = 10 GeV 2 with two different 
values of the cutoff z c on the light-cone momentum of the charmed quark (see 
text): the solid line is for z c = 0.5, the dot-dashed line for z c = 0.8. 
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